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ABSTRACT 

I investigate the Baryonic Tully-Fisher relation for a sample of galaxies with extended 21 cm rotation 
curves spanning the range 20 < Vf < 300 km s -1 . A variety of scalings of the stellar mass-to-light 
ratio T* are considered. For each prescription for T*, I give fits of the form Aid = <AVf. Presumably, 
the prescription that comes closest to the correct value will minimize the scatter in the relation. The 
fit with minimum scatter has A = 50 A4q km -4 s 4 and x = 4. This relation holds over five decades in 
mass. Galaxy color, stellar fraction, and are correlated with each other and with Aid, in the sense 
that more massive galaxies tend to be more evolved. There is a systematic dependence of the degree 
of maximality of disks on surface brightness. High surface brightness galaxies typically have T* ~ | 
of the maximum disk value, while low surface brightness galaxies typically attain 
Subject headings: dark matter — galaxies: kinematics and dynamics — galaxies: spiral 



j of this amount. 



1. INTRODUCTION 

The Tully-Fisher relation (Tully & Fisher 1977) has 
played an important role in establishing the extragalac- 
tic distance scale (e.g., Pierce & Tully 1988; Sakai et 
al. 2000). In this context, it has been treated as a sim- 
ple and convenient empirical relation between luminos- 
ity and line- width. The reason why it works is also im- 
portant, particularly to our understanding of galaxies as 
physical objects and how they formed (e.g., Eisentstein 
& Loeb 1996; McGaugh & de Blok 1998a; Courteau & 
Rix 1999; van den Bosch 2000; Navarro & Steinmetz 
2000a, b). Now that the issue of the distance scale is 
widely considered to be settled, we can hope to place an 
absolute scale on galaxy mass as well as luminosity. 

The physical basis of the Tully-Fisher relation is widely 
presumed to be a relation between a galaxy's total mass 
and rotation velocity (e.g., Freeman 1999). Luminosity 
is a proxy, being proportional to stellar mass, which in 
turn depends on the total mass. McGaugh et al. (2000) 
found that a more fundamental relationship between the 
baryonic mass and rotation velocity does indeed exist, 
provided that both stellar and gas mass are considered 
(Milgrom & Braun 1998). 

The relation resulting from the sum of stellar and gas 
mass is referred to as the Baryonic Tully-Fisher (BTF) 
relation. The BTF has also been investigated by Bell & 
de Jong (2001), Verheijen (2001), Gurovich et al. (2004), 
and Pfenniger & Revaz (2005). These efforts find broadly 
similar results, in that there is such a relation. How- 
ever, details of the relation differ. McGaugh et al. (2000) 
found a steep slope (x ~ 4) from a sample dominated by 
galaxies with H or /-band photometry and rotation ve- 
locities estimated from the line- width W20. Verheijen 
(2001) found much the same from if'-band photometry 
and flat rotation velocities measured from resolved 21 cm 
cubes. Both these authors assumed a constant value of 
the stellar mass-to-light ratio for all galaxies. Using Ver- 
heijen's data and a grid of stellar population models to 

Electronic address: ssm@astro.umd.edu 



refine the estimate of stellar mass, Bell & de Jong (2001) 
found a somewhat shallower slope (x « 3.5). Gurovich 
et al. (2004) find a break in the relation, with low mass 
galaxies following a steeper slope. 

There is an equal variety in the physical interpreta- 
tions. McGaugh et al. (2000) argue that the regularity of 
the BTF implies that, after the observed stars and gas are 
accounted for, no further comparably massive reservoirs 
of baryons are likely to exist in disk galaxies. Pfenniger 
& Revaz (2005) use the same data to argue the opposite 
case: the scatter is somewhat reduced if there are dark 
baryons weighing several times the observed gas mass. In 
the context of ACDM, one would expect still more dark 
baryons in order to match the universal baryon fraction 
(e.g., Mo & Mao 2004), though these need not be asso- 
ciated with the disk. 

The situation at present remains confused. The pur- 
pose of this paper is to provide the best empirical BTF 
relation possible with the currently available data. Much 
depends on the value of the stellar mass. I consider the 
effects on the BTF of varying the stellar mass over a 
broad range. Mass-to-light ratios are scaled by several 
recipes: as a fraction of maximum disk; with respect 
to stellar population synthesis models; and by the Mass 
Discrepancy — Acceleration (MDAcc) relation (McGaugh 
2004). This treats MOND (Milgrom 1983) as a purely 
phenomenological prescription. A grid of BTF fits are 
given that covers essentially any plausible choice of stel- 
lar mass-to-light ratio. The scatter of the relation varies 
with this choice, and an optimal choice that minimizes 
the scatter in the BTF relation is clear. 

2. THE DATA 

McGaugh et al. (2000) have already described the BTF 
over a large range of velocity and mass. This large dy- 
namic range is critical to constraining the slope of the 
relation, and also its absolute normalization since low 
mass galaxies are frequently gas rich and insensitive to 
assumptions about T*. Few other samples cover such a 
large dynamic range (e.g., Gurovich et al. 2004). 
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One thing that stands to be improved is the accuracy The use of Vf provides a considerable improvement in 

of the data. The data used here are from the sample of accuracy over line-width measurements (Verheijen 2001). 

Sanders & McGaugh (2002), as trimmed for accuracy by The result is a much cleaner BTF relation that is not af- 

McGaugh (2004). The reader is referred to these papers, flicted by possible corrections to line- widths for turbulent 

and references therein, for a further description of the motion. Turbulent corrections have the potential to af- 

sample. These galaxies all have extended 21 cm maps feet the slope of the BTF by systematically adjusting the 

from which the flat rotation velocity Vf is measured. line-width inferred rotation velocities of low mass galax- 
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ies. In addition to resolved atomic gas measurements, 
all galaxies have detailed surface photometry from which 
stellar masses can be estimated. 

The sample is identical to that in McGaugh (2004), 
with one exception. After re-examining the data, NGC 
2915 has been replaced by UGC 6818. Though the indi- 
vidual data points for NGC 2915 are quite precise, this 
galaxy does not have a well defined Vf, the quantity of in- 
terest for the BTF. The relative distance uncertainty for 
NGC 2915 is also uncomfortably large (Meurer, Mackie, 
& Carignan 1994; Karachentsev et al. 2003). In con- 
trast, UGC 6818 only barely missed the cut imposed in 
McGaugh (2004), and has more robust global measure- 
ments. Having done this exercise, I can see where further 
improvements could be made for individual galaxies, but 
these are generally very minor. These are all the galax- 
ies that are available with an obtainably high standard 
of accuracy. 

The data are given in Table 1. Column 1 gives the 
name of the galaxy. Column 2 gives Vf in kms -1 . These 
are the fit values given by Sanders & McGaugh (2002) 
which are the average of the outer points. As a test, I 
have remeasured Vf by eye for the sub-sample of Verhei- 
jen & Sancisi (2001). These agree to within a few kins" 1 
with the values given by Verheijen (2001) and with those 
tabulated here. Vf is easily and robustly measured, pro- 
vided only that the rotation curve is extended enough. 
Column 3 gives the stellar mass for the mass-to-light ra- 
tio from column 10, and column 4 gives the gas mass 
(both in units of 10 10 Mq). Column 5 gives the central 
surface brightness of the disk in B mag. arcsec -2 , and col- 
umn six the exponential disk scale length. The B-band 
is the only band-pass in common to all galaxies; consis- 
tent results are found in the subset with K'-h&nd data 
(Sanders & Verheijen 1998; Verheijen 2001; McGaugh 
2004). Column 7 gives the B — V color obtained from 
the original source given in Sanders & McGaugh (2002) 
if available, or from NED 1 if not. Colors measured with 
CCDs are given preference; if these are not available, 
RC3 values are used. B — V colors in square brackets are 
inferred from other measured colors through stellar pop- 
ulation models. It was most often the case that B — R 
had been measured instead of B — V . The colors are used 
to infer the mass-to-light ratio of the stellar population 
from said models. Several possible mass-to-light ratios 
are given in the last three columns: maximum disk in 
column 8, stellar population synthesis in column 9, and 
MDAcc in column 10 (McGaugh 2004). 

3. METHOD 
The BTF is expressed as 

M d = AVf, (1) 
where A is the normalization and x the slope. Note that 
since mass is used here rather than magnitudes, the slope 
m in traditional magnitude units would bem= — 2.5x. 
Vf is the measured rotation velocity in the flat part of the 
rotation curve, and A4d represents all measured 2 bary- 
onic mass. It includes both stars and gas. Most of these 

1 This research has made use of the NASA/IPAC Extragalactic 
Database (NED) which is operated by the Jet Propulsion Labo- 
ratory, California Institute of Technology, under contract with the 
National Aeronautics and Space Administration. 

2 In McGaugh ct al. (2000), this was referred to as the baryonic 
disk mass, hence the subscript. Both there and here the intent is to 



galaxies are disk dominated, but bulge mass is included 
where present (see Sanders 1996). The gas mass is that 
measured in neutral hydrogen, corrected for helium and 
metals. Other gas phases are presumed to be negligible. 
In effect, the molecular gas mass is subsumed into the 
stellar mass (see discussion in McGaugh 2004). By the 
same token, any net internal extinction is also subsumed 
by T*. This can not be a large effect, given the con- 
sistency between B and K' bands (Sanders & Verheijen 
1998). 

The baryonic mass of a galaxy is therefor 

M d = M*+M g . (2) 

The stellar mass is given by 

M* = T+L, (3) 

where T* is the mass-to-light ratio of the stellar popu- 
lation. Fig. 1 illustrates the various Tully-Fisher rela- 
tions that can be constructed from the data in Table 1: 
the luminosity- rotation velocity relation, L-Vf, this con- 
verted into stellar mass, M*-Vf, that with gas only (no 
stars), Mg-Vf, and the BTF, M d -V f . 

The choice of T* for each galaxy is critical. In Fig. 1 , 
the behavior in the A4*-Vf plane is markedly different 
from that in L-Vf. Use of the MDAcc mass-to-light ratio 
substantially reduces the scatter when mapping from L 
to M.+. In addition, the break at Vf « 90kms~ 1 noted 
by McGaugh et al. (2000) becomes apparent in A4*-Vf, 
though it is not visible in L-Vf. Constant mass-to- light 
ratios were used in McGaugh et al. (2000), but the break 
was apparent there because of the large number of very 
low mass, gas dominated dwarfs. The data quality re- 
strictions imposed here exclude those objects. Simply 
applying a constant Y* to the galaxies here would pre- 
clude the discovery of the break, since this would just be 
a shift in the scale of L-Vf. Yet a specific, well-defined 
method for determining Y* recovers the break without 
input about its existence. 

I explore three choices for the stellar mass-to-light ra- 
tio: maximum disk, T max ; stellar population synthesis 
models, Y pop ; and that from the MDAcc, Y acc . For each 
of these choices, I construct a grid of Y* scaled from each 
by a constant factor: T, V, or Q: 

Y* = TT max (4) 
Y* = PT pop (5) 
T* = QY acc . (6) 

This provides a prescription for Y* that is specified by 
method and scaling factor. Fig. 2 illustrates the BTF 
relations stemming from various choices. A few details 
concerning each method are worth noting here. 

3.1. Scaling by Maximum Disk 

Maximum disk is the highest mass-to-light ratio con- 
sistent with but not exceeding the rotation curve data. 
The inner shape of rotation curves are often quite con- 
sistent with the shape predicted by the observed baryons 
(van Albada & Sancisi 1986; Selwood 1999; Palunas & 
Williams 2000), leading some to argue that disks must 
be nearly maximal. I consider the range < T < 1. 

represent all known, observed baryons, regardless of their location 
in the galaxy (bulge or disk). 
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Fig. 1. — Four versions of the Tully-Fisher relation. The top left panel shows the B-band luminosity as a function of the flat rotation 
velocity. The top right panel plots the gas mass instead of luminosity. The bottom left panel plots stellar mass, J\4* = T*L for the MDAcc 
mass-to-light ratios. The bottom right panel plots the Baryonic Tully-Fisher relation, with Md = M* + Mg- The solid line is a fit to the 
data, M d = 50V^. This is drawn as a dashed line in the other panels for comparison. In the top left panel, the line is drawn for the mean 
value <T*) = 1.7 Mq/Lq. 



I adopt the maximum disk value given by the origi- 
nal source for the data for each galaxy (see Sanders & 
McGaugh 2002). Caution should be exercised in inter- 
preting T max - Some authors leave room for a dark halo, 
so that even "maximum" disk may provide only ~ 84% 
of the observed velocity at the peak of the disk contri- 
bution (Sackett 1997). Others leave zero room, fitting 
disk-only models as far out as possible (e.g., Palunas & 
Williams 2000). This is a fairly subtle distinction in ro- 
tation curve decompositions, but does make a noticeable 
difference in Y*. For example, the average 7-band mass- 
to-light ratio of the Palunas & Williams (2000) sample is 
{^max) = 2AMq/ 'Lq. Since M cx V 2 , if we scaled this 
down from a 100% to 84% contribution, the mean would 
be {T max ) = 1.7 AAq/Lq (see also Barnes, Sellwood, & 
Kosowsky 2004). 

In general, there is no uniform definition of maximum 
disk. The data originate from a wide variety of sources, 
so there is no guarantee as to how maximal maximum 
disk is. Moreover, things can go both ways. For low 
surface brightness disks, very small differences in the ro- 
tation curve can lead to large changes in the inferred 
value of T max (Swaters et al. 2000; McGaugh, Rubin, & 
de Blok 2001). For these galaxies and some of the more 
extreme dwarf galaxies, I have re-assessed the value of 
maximum disk. In spite of these caveats, the original 
value appears sensible for most of the galaxies in Table 
1. 



There are a few exceptions. For the Ursa Major data 
(Verheijen & Sancisi 2001), which comprises a substan- 
tial plurality of the data in Table 1, a "soft" limit was 
imposed on the value of T max so that Vj of the fit- 
ted pseudo-isothermal halo did not diverge to absurdly 
large values (Verheijen 1997). This can easily happen 
for galaxies that are well described by maximum disk 
and for which the data are not very extended in radius. 
There is little clear need for a halo in such cases (Palu- 
nas & Williams 2000), so Vf is poorly constrained and 
tends towards large values since only the rising portion 
of the contribution of the dark matter is seen. Unfor- 
tunately, imposing this "soft" constraint can sometimes 
lead to T max which is very sub-maximal. This is a subtle 
point which only became apparent because T max < T acc 
for four of the galaxies in Table 1 (Sanders & Verheijen 
1998). This is a mathematical impossibility, so I have set 
^ max = Y acc in these cases. 

3.2. Scaling from Stellar Population Synthesis Models 

Stellar population models have advanced to the point 
where they give plausible estimates of the mass-to-light 
ratio, even for the composite stellar populations of spiral 
galaxies. They are not yet perfect of course, but do pro- 
vide a decent choice for estimating Y* (Bell & de Jong 
2001; Portinari et al. 2004). Here I employ the models 
of Bell et al. (2003) to estimate the stellar mass-to-light 
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Fig. 2. — The BTF for various choices of stellar mass-to-light ratio. The left column of panels shows scalings relative to maximum disk: 
r = 1, 0.5, and 0.25 from top to bottom. The middle column shows scalings relative to the population synthesis models of Bell ct al. (2003) 
V = 2, 1, and 0.5 from top to bottom. Similarly, the right column shows scalings relative to the mass-to-light ratio from the MDAcc 
Q = 2, 1, 0.5 from top to bottom. (Note that r = V = Q = are all equivalent to the gas-only panel in Figure 1.) The mass-to-light ratio 
is not allowed to exceed the maximum disk value. Galaxies arc plotted as open symbols with their mass-to-light ratios set to the maximum 
disk value if the value specified by V or Q would have exceeded maximum disk. Half of the sample has reached this point by Q = 2. The 
fit to the Q = 1 case is shown as a dashed line in all panels for comparison. 



ratio from the observed color: 

log T pop = 1.737(5 - V) - 0.942 (7) 

(from their Table 7). 

The B — V color is given precedence in estimating T pop . 
When B — V is not available, whatever color is available 
is used. The most common substitute is B — R, which 
is expected to be nearly as well correlated with T* as 
B — V (Bell & de Jong 2001). The same model (from 
Table 7 of Bell et al. 2003) is used to estimate B - V 
(the bracketed colors in Table 1), but T pop is based on 
the observed color. Credible color information could not 
be located for UGC 2259 and DDO 170, the only pieces 
missing from Table 1. 

3.3. Scaling from the Mass-Discrepancy — Acceleration 

Relation 

McGaugh (2004) used the detailed shapes of the ro- 
tation curves of the galaxies in Table 1 to show that 
there is an empirical relation between acceleration and 
the amplitude of the mass discrepancy (essentially the 
ratio of dark to baryonic mass). This relation holds at 



every point along a resolved rotation curve for any non- 
zero choice of mass-to-light ratio. The scatter about this 
relation depends on this choice; T acc is determined by 
minimizing the scatter with respect to the mean local re- 
lation. It is interesting to see here how this T* estimator 
fares with the global BTF. 

The MDAcc is a purely empirical relation. It is mathe- 
matically equivalent to MOND: the T acc are the same as 
the MOND best fit values (Begeman, Broeils, & Sanders 
1991; Sanders 1996; Sanders & Verheijen 1998; de Blok 
& McGaugh 1998). Here we must make the distinction 
between MOND as a fundamental theory (with its asso- 
ciated difficulties), and as a successful recipe for fitting 
rotation curves. Only the latter is required. Indeed, the 
MDAcc is the local analog of the BTF, and had MOND 
never been invented, we would perhaps already have rec- 
ognized the MDAcc purely as an empirical relation (San- 
cisi 2003). Just as the Tully-Fisher relation can be used 
empirically to estimate distances without understanding 
its physical basis, so too can the MDAcc be utilized to es- 
timate stellar mass-to-light ratios without prejudice con- 
cerning its theoretical basis. 
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4. RESULTS 

Equations (4-6) are used to estimate the stellar mass 
for the various scalings. For each scaling, a grid of 10 
choices of the scaling constant are made: T = 0.1 to 1.0 
in steps of 0.1 relative to maximum disk, and for V and 
Q values ranging from 0.2 to 2.0 in steps of 0.2. The 
maximum disk scaling obviously can not exceed unity, 
while for the other scalings there is no reason not to 
consider values larger than one. For example, V > 1 
would simply imply an IMF heavier than assumed in 
the nominal population model which has been adopted. 
However, these values should not exceed maximum disk, 
and arc not allowed to do so. If the choice of V or Q 
exceeds the maximum disk value for a particular galaxy, 
the maximum disk value is used instead. Half of the 
sample has saturated at maximum by Q = 2, so larger 
values are not considered. 

The BTF (logarithm of equation 1) is fit 3 to the data 
in Table 1 for each set of choices for T*. Every galaxy 
carries equal weight in the fits. Since the data have been 
selected to be of high quality, the residual uncertainties 
are likely to be dominated by systematic effects (such as 
the precise distance to each galaxy) rather than factors 
internal to the data, though these obviously matter as 
well. In any case, the choice of mass-to-light ratio com- 
pletely dominates the results. 

The results of fitting the BTF are given in Table 2. 
The intercept log .4 and slope x are recorded, together 
with the formal uncertainties in each (a a and a x ). Also 
given is the scatter of the data about each relation, um . 
For reference, the input data, in the form of the normal 
luminosity-based Tully-Fisher relation is given. It has a 
scatter ul — 0.24, which is equivalent to 0.6 mag. (base 
ten logarithms are used throughout). Note also that the 
limit T = 'P=Q = 0is equivalent to the gas-only rela- 
tion, so is given only once. 

The results summarized in Table 2 are illustrated in 
Figs. 3 and 4. As the assumed mass-to-light ratio in- 
creases, the slope x gradually increases while the zero 
point A decreases to compensate. The lowest mass galax- 
ies in Table 1 are dominated by gaseous rather than stel- 
lar mass, so the BTF tends to pivot about them. This 
can be seen by eye in Fig. 2, where the low end of the 
relation hardly budges while the more massive, star dom- 
inated galaxies move up and down with Y*. These low 
mass, gas-rich galaxies provide a critical anchor point for 
the absolute calibration of the BTF since their location 
in this diagram is insensitive to the choice of T*. 

The scatter cjm about each fit is shown in Fig. 4. 
The maximum disk scaling has scatter comparable to 
the scatter in the input luminosities. The use of color 
information with population synthesis models provides 
a small increment of improvement in the scatter, as it 
should if the models succeed in improving the estimate 
of Y* over a constant value for all objects. As stellar 
mass is rendered unimportant for very small T and V 
the scatter starts to increase. The limit of zero stellar 
mass, with gas only, obviously makes for an inadequate 
BTF, as one would expect. 

4.1. The Optimal BTF 

3 These are 'direct' fits. The scatter is small enough that forward 
and reverse fits are not distinguishable except in the limit Y* — > 0. 



TABLE 2 
BTF Statistics 





log .4 




X 






T T, 


Luminosity (Normal TF) 






1.0 


2.52 


0.37 


3.48 


0.17 


0.239 


A/f 

JVlg 


Gas Only 










0.0 


4.67 


0.48 


2.28 


0.48 


0.313 


r 


Maximum 


Disk Scaling 








0.1 


3.24 


0.37 


3.09 


0.37 


0.244 


0.2 


2.76 


0.35 


3.38 


0.35 


0.230 


0.3 


2.50 


0.35 


3.55 


0.35 


0.224 


0.4 


2.34 


0.34 


3.67 


0.34 


0.222 


0.5 


2.23 


0.34 


3.75 


0.34 


0.220 


0.6 


2.15 


0.34 


3.81 


0.34 


0.220 


0.7 


2.09 


0.34 


3.87 


0.34 


0.220 


0.8 


2.05 


0.34 


3.91 


0.34 


0.220 


0.9 


2.01 


0.34 


3.94 


0.34 


0.220 


1.0 


1.98 


0.34 


3.97 


0.34 


0.221 


V 


Population Synthesis Scaling 






0.2 


3.81 


0.33 


2.81 


0.15 


0.204 


0.4 


3.44 


0.30 


3.04 


0.14 


0.186 


0.6 


3.23 


0.29 


3.19 


0.14 


0.183 


0.8 


3.08 


0.29 


3.29 


0.13 


0.183 


1.0 


2.98 


0.29 


3.37 


0.13 


0.186 


1.2 


2.87 


0.28 


3.44 


0.13 


0.188 


1.4 


2.75 


0.27 


3.52 


0.13 


0.191 


1.6 


2.65 


0.27 


3.58 


0.13 


0.194 


1.8 


2.54 


0.26 


3.64 


0.12 


0.197 


2.0 


2.46 


0.26 


3.69 


0.12 


0.200 


Q 


MDAcc Scaling 








0.2 


3.05 


0.25 


3.18 


0.25 


0.165 


0.4 


2.46 


0.20 


3.53 


0.20 


0.129 


0.6 


2.11 


0.17 


3.74 


0.17 


0.112 


0.8 


1.87 


0.16 


3.89 


0.16 


0.103 


1.0 


1.70 


0.15 


4.00 


0.15 


0.098 


1.2 


1.60 


0.15 


4.07 


0.15 


0.100 


1.4 


1.58 


0.16 


4.10 


0.16 


0.107 


1.6 


1.60 


0.18 


4.10 


0.18 


0.118 


1.8 


1.61 


0.20 


4.10 


0.20 


0.129 


2.0 


1.64 


0.22 


4.10 


0.22 


0.140 



Irrespective of the physics underlying the BTF, we 
expect that the prescription for T* that comes closest 
to the correct value will minimize the scatter about it. 
There may be some intrinsic scatter, but if we can get T* 
right, there should be no scatter left due to it. The T acc 
prescription gives less scatter than either T max or T pop . 
There is clear variation in the scatter with Q, with a well 
defined minimum at Q = 1. The same mass-to-light ratio 
chosen to minimize the residuals from the local MDAcc 
also minimizes the scatter in the global BTF. This was 
already apparent in Fig. 6 of McGaugh (2004), and is 
quantitatively confirmed here. It is rare in extragalac- 
tic astronomy that we find a correlation as strong as the 
BTF (K = 0.99 for Q = 1). Indeed, the BTF is so 
sharply defined that it has a kurtosis of 1.5. 

It is difficult to avoid the conclusion that the choice 
Q = 1 does effectively give the correct Y*. This provides 
an absolute calibration of the BTF: 

M d = 50V/ (8) 

with Aid in A4q and Vf in kms -1 . Indeed, the scatter 
in this BTF relation is so small that it could be used 
to estimate Y* with nearly as great accuracy as the full 
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r p q 

Fig. 3. — The slope (top panels) and zero point (bottom panels) of the BTF for various scalings of the stellar mass-to-light ratio: as a 
fraction of maximum disk (left panels), relative to population synthesis models (middle panels), and relative to the MDAcc mass-to-light 
ratios (right panels). For each choice of Y*, the data have been fit as a straight line to the logarithm of M d = -AVJ (Table 2). The thin 
lines show the formal uncertainties on the parameters log A and x of the fit. For all scalings, the slope of the mass — circular velocity BTF 
only becomes as shallow as x = 3 for implausibly small mass-to-light ratios. 




O 1 1 

0.5 1 1.5 2 
2T, P, Q 

Fig. 4. — The scatter aj^ m the BTF for the various scalings of 
mass-to-light ratio. Dashed line: maximum disk (plotted against 
2r to fill the same range as the other choices). Dotted line: popu- 
lation synthesis. Solid line: mass-to-light ratios from the MDAcc. 
The latter show a clear minimum at Q = 1. The other scalings 
have no clearly defined minima, and larger scatter for all choices of 
r or V. The ordinary Tully-Fishcr relation of the input data has 
a scatter of ox = 0.24 (dash-dot line). 



MDAcc, but with considerably less information. To ap- 
ply the latter, we require a well resolved rotation curve 
and Hi surface density map, and detailed surface pho- 
tometry. The BTF requires only global quantities: L, 



A4 g from a single dish 21 cm observation, and an esti- 
mate of Vf (preferably from a rotation curve which is suf- 
ficiently resolved to perceive the flat part, though 
would be an adequate substitute, with some penalty in 
accuracy.) 

This BTF is consistent with the relation first reported 
by McGaugh et al. (2000). The slope is identical, but 
the normalization is somewhat larger. This change is 
not particularly significant (lcr). It is largely due to the 
improvement in data quality and the use of Vf rather 
than jW 2 a- The line- width W 2 q is systematically larger 
than 2Vj, which is more closely approximated by W 50 
(Broeils 1992; Verheijen 1997). Thus we should expect 
some increase in A simply from the change of circular 
velocity measures. 

Indeed, it is reassuring that the relation found here 
is so closely consistent with the previous version. Mc- 
Gaugh et al. (2000) chose Y* by a different method, tak- 
ing a constant value for all galaxies motivated by popu- 
lation models. More importantly, the data are largely 
independent: the BTF of McGaugh et al. (2000) in- 
cluded 110 galaxies from Bothun et al. (1985), 14 from 
Matthews, van Driel, & Gallagher (1998), and 65 from 
Eder & Schombert (2000) that are completely distinct 
from and independent of the data in Table 1 . The agree- 
ment of such diverse and independent data sets provides 
strong confirmation of the basic empirical result. 

4.2. Correlations Among Evolutionary Parameters 
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TABLE 3 
Correlation Coefficients 





M d 


/* 


T* 


B-V 


Vf 
M d 

f* 
T* 


0.99 


0.78 
0.76 


0.77 
0.76 
0.84 


0.53 
0.53 
0.54 
0.64 



The BTF is by far the strongest correlation present in 
the data, but it is not the only one. There are a few 
others worth noting. Some of these are shown in Fig. 5. 
Many others are not shown to avoid redundancy: any- 
thing that is correlated with disk mass is also correlated 
with Vf through the BTF. Correlation coefficients for the 
quantities appearing in Fig. 5 are given in Table 3. 

One interesting quantity that can be constructed from 
the data in Table 1 is the fraction of the total observed 
baryonic mass in the form of stars: /* = M.+j M. d . The 
variation of /* with disk mass basically shows the turn- 
away of the stellar mass TF from the BTF as one goes 
from star to gas domination. There is lots of scatter in 
the f±-Md diagram which is not reflected in the BTF. 
This reiterates that what matters is mass, not the form 
it is in. 

Together with the mass-to-light ratio and color, /* is an 
indicator of the evolutionary state of a galaxy (McGaugh 
& de Blok 1997, 1998a; Schombert, McGaugh, & Eder 
2001). As evolution proceeds and a galaxy converts its 
initial gas into stars, the stellar fraction increases. As 
old stars accumulate, the mean color reddens, and the 
mass-to-light ratio increases. One therefor expects these 
quantities to be related. 

The evolutionary quantities Y*, and color do indeed 
correlate as expected. Not only does the mass-to-light ra- 
tio increase as the stellar population reddens, so too does 
Y* climb as /* increases. Indeed, the Y*-/* correlation 
is the next strongest in Table 3 after the BTF. 

In this sample, the evolutionary quantities are corre- 
lated with baryonic mass. They are also correlated with 
Vf, as shown previously by Sanders (1996). This must 
follow, given the BTF. Indeed, figures constructed with 
Vf instead of Md look so similar that they are difficult 
to distinguish. 

One curious feature of the correlations with disk mass 
is that all of the evolutionary indicators suggest that 
more massive disks are typically more evolved. This is 
opposite what one might naively expect from a hierarchi- 
cal galaxy formation picture. In such a scenario, small 
galaxies form first so should have evolved the most. That 
the opposite is true seems more consistent with mono- 
lithic or even top-down galaxy formation. Similar results 
have been found recently at high redshift (e.g., Juneau 
et al. 2005; Trcu ct al. 2005), where the tendency for 
massive galaxies to be the most evolved is called "down- 
sizing." If galaxies do form hierarchically, the observed 
trend suggests that the mechanisms which regulate a 
galaxy's post-formation evolution dominate over the for- 
mation epoch in determining its present evolutionary sta- 
tus. This would appear to be a continuous function of 
mass, given the continuity in Fig. 5. It does not appear to 



be as simple as a single epoch of reionization suppressing 
galaxy formation at a characteristic mass scale. 

One word of caution is that while these data span a 
broad range of galaxy properties, they do not constitute 
a complete volume limited sample. The details of these 
correlations may well change as other data are added. 
Nevertheless, it seems unlikely that the trend apparent 
in this figure could be reversed: massive galaxies still in 
the early stages of evolution seem to be very rare. 

It is extremely difficult to obtain a sample extending 
to very low mass which is complete in any meaningful 
sense. Since an Hi map is a prerequisite for membership 
in this sample, it seems likely that there are some low 
mass, high /* galaxies which are not represented here. 
So while there are some clear correlations with disk mass 
in this sample, there is no guarantee that precisely these 
correlations will hold for all galaxies. There is more hope 
that correlations between the evolutionary parameters 
/*, Y*, and B — V will hold, if only because they make 
so much sense. 

These cautions do not apply to the BTF. There is no 
hint of deviation from it: the residuals are small and un- 
corrected with any measured parameter. There is no 
reason to suspect that any rotating galaxy deviates from 
the BTF — even those that should deviate do not (Mc- 
Gaugh & de Blok 1998a; Verheijen 2001). The BTF 
appears to be a fundamental relation. 

4.3. Consistency with Stellar Population Models 

The remarkable consistency of the MDAcc mass-to- 
light ratios with stellar population synthesis models has 
been noted previously (Sanders 1996; McGaugh & de 
Blok 1998b; Sanders & McGaugh 2002; McGaugh 2004). 
I have added more color information to the data in Table 
1 than has previously been available, and this point re- 
mains true. Indeed, there are now enough data that we 
can perform a fit in the same manner as done for popu- 
lation synthesis models. Following the format of Bell & 
de Jong (2001), the data are fit to a relation of the form 

logT* =a-b(B-V). (9) 

Portinari et al. (2004) note that while their £?-band 
model Y* are linearly correlated for B — V > 0.55, they 
show a break at this point. Y* turns down to lower val- 
ues for bluer colors than predicted by extrapolation of 
the line fit to redder colors. I therefor consider two fits 
to the data: one covering data of all colors, and the other 
restricted to B - V > 0.55 

The relations from the stellar population synthesis 
models of Bell et al. (2003) and Portinari et al. (2004) 
are compared with the fits to the data in Table 4. The 
relation of Bell et al. (2003) is that used to estimate 
Yp p (equation 7). Note the close agreement between 
the two models: these are virtually indistinguishable in 
Fig. 5. The zero points a of the models are dominated 
by the choice of IMF, which is fortuitously close (Kroupa 
& Weidner 2003). 

There is close agreement between the parameters a and 
b fit to the data and those predicted by the models. The 
MDAcc mass-to-light ratios, determined by dynamical 
methods completely independent of the population mod- 
els, show precisely the same trends. A slight offset in 
normalization is apparent, but its formal significance is 
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Fig. 5. — Some of the stronger correlations in the current sample beyond the BTF. The stellar mass fraction /* is correlated with disk 
mass and Y* (left panels). The mass-to-light ratio is also correlated with disk mass and color (right panels). In the T* — (B — V) diagram 
at bottom right, a fit to the data is shown as a solid line. The fit shown is to data with B — V > 0.55, where one expects a break in the 
T* — (B — V) relation (Portinari et al. 2004). The dashed line has the slope expected for B — V < 0.55 in their models, with normalization 
chosen to match the fit to the data at B — V = 0.55. For comparison, the dotted line shows the model mass-to-light ratio — color relation 
of Bell ct al. (2003) which has been used to define V = 1. 



TABLE 4 

Mass-to-Light Ratio-Color Relations 



Source 


a 


b 


Note 


Bell et al. 


-0.942 


1.737 


Scaled Salpeter IMF 


Portinari et al. 


-0.925 


1.69 


Kroupa IMF 


Fit to the Data 


-0.90 


1.83 


B-V > 0.55 


Fit to the Data 


-1.22 


2.31 


All colors 



Note. — Relations of the form 
logT* =a + b(B-V). 



low. It is tempting to interpret this offset as the molecu- 
lar gas which has been subsumed into T acc . If the mass 
of gas in the molecular phase is typically 10% to 20% of 
that in stars, the offset between models and data would 
be reconciled. 

The break at/? — y«0.55in the models of Portinari 
et al. (2004) is apparent in the data. Fig. 5 shows (as 
the dashed line) the model slope for bluer colors (from 
their Fig. 28), normalized to match the fit to the data 
at the break point. This is entirely consistent with the 
downward trend in the data. 

The scatter in Y* is also as expected, being larger in B 
than in K 1 (Sanders & McGaugh 2002; McGaugh 2004). 
This is true also above and below the break point at 



B — V = 0.55. Below this point, one expects a tremen- 
dous amount of scatter from the rapid evolution of young 
populations. Here the scatter is enormous: ux, = 0.52. 
Above this color, one expects variation in to settle 
down, as the effects of individual star formation events 
are moderated by the accumulation of mass from previ- 
ous generations. For these redder colors, the scatter is 
<7T 4 = 0.20. 

In sum, the Q = 1 mass-to-light ratios are optimal not 
only in terms of minimizing the scatter in the BTF and 
the MDAcc (from which they come), but also in terms 
of our expectations for stellar populations. Indeed, it is 
hard to imagine better agreement with independent pop- 
ulation models to which no fit has been made. Moreover, 
the scatter in the dynamical relations is so small for these 
high quality data that Y* estimated from the BTF itself 
(equation 7) are nearly indistinguishable from T acc from 
the MDAcc. Either empirical method can be employed 
with unprecedented precision. 

4.4. Test by Extrapolation 

The value of the slope x of the BTF is somewhat 
controversial, being of considerable physical importance. 
The nominal expectation of CDM is that x — 3 (Navarro 
& Steinmetz 2000a, b), and Courteau et al. (2003) argue 
that this is consistent with their /-band data. The I- 
band is a good but not perfect indicator of mass, and 
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the HST calibrated distances of Sakai et al. (2000) give 
an /-band slope of 4. This emphasizes the need for a 
good estimate of Md and not just L. 

McGaugh et al. (2000) found x = 4 from data and 
Y* estimates independent of those used here. The key 
aspect of that study was the large number of very gas 
rich, low rotation velocity galaxies which tied down the 
low mass end of the BTF. Verheijen (2001) obtained the 
same result. Bell & de Jong (2001) generated stellar pop- 
ulation models in an attempt to improve the estimate of 
T*. Applied to the data of Verheijen, they found a BTF 
with x — 3.5. I have, in effect, repeated this analysis with 
more data and updated models, and find very much the 
same result: for V = 1, x = 3.4 (Table 2). 

It is interesting that while the models and data are in 
good agreement (§4.2), use of T pop gives somewhat shal- 
lower x than T acc or T max (Fig. 3). This is attributable 
in part to the slight difference in normalization (a in 
Table 4). This can not be the entire reason for the dif- 
ference, as increasing V to 2 only increases x to 3.7. The 
shortcoming of T pop estimated from a single color is that 
galaxies become beads on a string in the Y+-(B — V) 
diagram: all the data collapse to fall on the dotted line 
in Fig 5. The inability to estimate deviations from the 
mean relation and thereby include a realistic estimate of 
the scatter seems to result in a bias 4 in the determina- 
tion of the slope. Relative to the population predicted 
Tfp p, T acc is skewed to high values for red colors and 
to low values for blue colors. One expects such a skew 
relative to a single straight line fit (Portinari et al. 2004) ; 
it causes a systematic difference in the slope in spite of 
the close agreement seen in Fig. 5 and Table 4. 

The critical issue for constraining the slope is the dy- 
namic range of the data. As a further test of the slope 
of the BTF, I have sought out galaxies that extend the 
relation to lower rotation velocities. The slowest rotator 
in Table 1 has Vf = 54kms~ 1 . I have searched the litera- 
ture for galaxies with slower rotation velocities that are of 
adequate quality to make a useful comparison. There are 
many slow rotators already in the sample of McGaugh et 
al. (2000), but there the rotation velocity estimate was 
based on a line- width (Eder & Schombert 2000). Here I 
require that there be a resolved measurement adequate 
for estimating Vf. 

Table 5 contains the objects found that met these cri- 
teria. Many of these are due to the recent excellent work 
of Begum & Changular (2003, 2004a,b) which extends 
down to objects approaching the globular cluster mass 
scale. In Table 5 I give the data necessary for the BTF 
from the data given in each reference cited there. I have 
estimated Vf myself from the published data. In ad- 
dition to the best estimate, I give a generous range of 
uncertainty. In many of these cases, there is a substan- 
tial correction for asymptotic drift. The lower limit of Vf 
takes the lower end of the published errors ignoring the 
asymptotic drift correction. The upper limit takes the 
upper end of the experimental errors and includes the 
correction. This procedure results in a very broad error 
estimate: the uncertainties listed in Table 5 are much 

4 One could, perhaps, do better with more colors. The incorpo- 
ration of such information had an analogous effect on the slope of 
the if'-band Y*-color relation, which became flatter between Bell 
& dc Jong (2001) and Bell et al. (2003). 



TABLE 5 
Extreme Dwarf Galaxy Data 



Galaxy 


Vf 


M* Mg 


Ref. 




(kms- 1 ) 


(10 6 Mq) 





ESO215-G?009 


51 I* 


23 


714 


1 


UGC 11583 a 


48 ±\ 


119 


36 


2,3 


NGC 3741 


44 +_t 


25 


224 


4 


WLM 


38 t\ 


31 


65 


5 


KK98 251 


36 t\ 


12 


98 


3 


GR 8 


25 +_\ 


5 


14 


6 


Cam B 


2oi!° 


3.5 


6.6 


7 


DDO 210 


17 ±1 


0.9 


3.6 


8 



References. — 1. Warren, Jerjen, & Koribalski (2004). 
2. McGaugh, Rubin, & de Blok (2001). 3. Begum & Chen- 
galur (2004a). 4. Begum, Chengalur, & Karachentsev (2005). 
5. Jackson et al. (2004). 6. Begum & Chengalur (2003). 7. Be- 
gum, Chengalur, & Hopp (2003). 8. Begum, & Chengalur 
(2004b). 

a UGC 11583 = KK98 250. 



larger than la. 

For the stellar and gas mass, I take the value given by 
the original authors. has often been estimated using 
the models of Bell & de Jong (2001), so these correspond 
roughly to V = 1 estimates. Rather than an uncertainty, 
I consider the full range of possible stellar masses, from 
zero to maximum disk. These are the vertical lines in 
Fig. 6. 

The extrapolation of the BTF fit to the data in Table 
1 does an excellent job of predicting the data for these 
lower mass galaxies. The objects in Table 5 follow the 
slope x — 4 down to unprecedented low velocity. The 
BTF remains valid over five decades in mass. 

Fig. 6 also illustrates why the slope has been difficult to 
constrain. The lower limit of other studies (~ 10 9 Mq, 
V f w 70 kms" 1 : e.g., Bell & de Jong 2001; Courteau 
et al. 2003) is shown in the left hand panel. The data 
considered there extend over only a small fraction of the 
range studied here (down to 5 x 10 6 A4q with the extreme 
dwarfs in Table 5). This lack of dynamic range probably 
dominates all other factors (such as the mass determina- 
tion method or the velocity estimator used) in constrain- 
ing the slope. Indeed, if we were to truncate the data 
in Table 1 at the same level, we would fail to perceive 
the break-point in the stellar mass Tully-Fisher relation. 
Such a sample would be dominated by star-dominated 
galaxies, and fail to provide the constraint on Y* which 
follows when gas-dominated objects are included. 

The steep slope of the BTF should come as no sur- 
prise, as it is completely consistent with the results of 
McGaugh et al. (2000). That study made use many 
low mass galaxies, 19 of which have 5W20 < 70kms~ 1 . 
Those objects are completely independent of the galax- 
ies discussed here. The new, more accurate data in Table 
5 simply return the same result with less scatter. Other 
workers investigating low mass galaxies have also inferred 
the need for a steep slope (e.g., Gurovich et al. 2004; 
Pizagno et al. 2004). 

Formally, a slope as shallow as x = 3 deviates from the 
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Fig. 6. — The stellar mass (left) and baryonic (right) Tully-Fisher relations, including the data for the extreme dwarf galaxies listed in 
Table 5. The horizontal lines through these objects are the maximum plausible range for Vf. these are much larger than la error bars. 
The vertical lines show the full range of possible stellar masses, from zero to maximum disk. The extrapolation of the BTF fit to the more 
massive galaxies from Table 1 is in good agreement with these extreme dwarfs. The importance of this check is illustrated by the thin lines 
inset in the left panel. These show the limits of samples that suggest shallower slopes (e.g., Bell & de Jong 2001; Courteau et al. 2003). 



optimal BTF by la. This can be made less by changing 
Y*, but only at the price of degrading the correlation 
and the the many consistency checks on T*. In order to 
recover a slope as shallow as x = 3, one requires V = 0.36 
or r < 0.1 (Table 2). Such absurdly sub-maximal disks 
would fall considerably short of the mass which is directly 
observed in stars locally. 

Consideration of the extreme dwarfs renders it even 
more difficult to reconcile a shallow slope with the data. 
While it is possible, at least in principle, to move massive 
galaxies down in mass by reducing their mass-to-light 
ratios in order to accomodate a shallow slope, it is not 
possible to move the extreme dwarfs very far up. Even 
taking maximum disk in those cases makes little differ- 
ence to the slope, and causes the curious situation that 
the IMF in these low mass galaxies must be systemati- 
cally heavier than that in giant galaxies. It is thus very 
difficult to reconcile a shallow [x = 3) slope for the BTF 
with the data. 

4.5. The Maximality of Disks 

One application of the result here is to quantify the de- 
gree to which galaxy disks are maximal. There is consid- 
erable debate as to whether high surface brightness disks 
are maximal (e.g., Sellwood 1999; Courteau & Rix 1999). 
There would seem little doubt that low surface brightness 
disks are dark matter dominated (de Blok & McGaugh 
1997, 2001), but an argument for maximal disks can be 
made even in these objects (Fuchs 2003). It is therefor 
of considerable interest to investigate how maximal disks 
are, and how disk maximality varies with disk properties. 

For the MDAcc mass-to-light ratios favored here, the 
fraction of maximum disk in each case is 

r* = -^l. (io) 

1 max 

This is plotted against disk mass and surface density 



in Fig. 7. There is only a weak correlation of with 
disk mass (1Z = 0.45) which depends heavily on rather 
few points at low mass (cf. Persic & Salucci 1988). Dy- 
namical arguments stemming from the adherence of low 
surface brightness galaxies to the Tully-Fisher relation 
(Zwaan et al. 1995; Sprayberry et al. 1995; Hoffman et 
al. 1996) suggest that disk maximality should cor- 
relate with surface brightness (Tully & Vcrheijen 1997; 
McGaugh & de Blok 1998a; Zavala et al. 2003). We can 
improve on this by using the mass-to-light ratios T acc to 
convert the observed central surface brightness into the 
central surface mass density of stars: 

log S - logT acc + 0.4(27.05 - fi ). (11) 

As anticipated, there is a good correlation between T* 
and S (1Z = 0.74). A fit to the data in Fig. 7 gives 

log r* = -0.98 + 0.3 log S . (12) 

In terms of the more directly observable central surface 
brightness, this translates to logT* = 3.13 — 0.16/io- 
There is considerably greater scatter about this latter 
relation. 

Irrespective of how we frame the relation, or what 
mass-to-light ratio prescription we prefer, it seems in- 
evitable that the disk contribution must decline system- 
atically as surface density declines. Low surface bright- 
ness disks are inevitably dark matter dominated. In 
contrast, high surface density disks contribute a non- 
negligible fraction of the total mass at small radii for 
plausible T*. Remarkably, this leaves no residual signa- 
ture in the Tully-Fisher relation (McGaugh & de Blok 
1998a; Courteau & Rix 1999) in spite of the generally 
modest radius at which rotation curves achieve Vf. 

Statistics of these data, divided into quartiles by S , 
are given in Table 6. The typical = 0.78 in the high- 
est surface density quartile. can not exceed unity, 
and is projected to saturate at fi ~ 19.5 mag. arcsec -2 . 
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Fig. 7. — The maximality of disks (T* for Q = 1) as a function of disk mass (left) and the central stellar surface mass density (right). 
While the correlation with mass is poor, there is a clear correlation of maximality with disk surface density. High surface density disks 
tend to be nearly maximal, while lower surface density disks are systematically sub- maximal. The line is a fit to the data (equation 12). 



TABLE 6 
Maximality of Disks 



Quartile 


N 


(mag. arcsec 2 ) 


(So) 
(M oP c- 2 ) 


<r*> 


1 


15 


23.31 


21 


0.25 


2 


15 


22.31 


99 


0.48 


3 


15 


21.52 


295 


0.74 


4 


15 


20.41 


964 


0.78 



Note. — The biweight location of the the disk central 
surface brightness, central stellar mass surface density, and 
the degree of maximality of disks are given for each quartile 
of the sample. 



This is comparable to the highest surface brightness disks 
that exist (Marshall 2004). In the lowest surface density 
quartile, the typical fraction of maximum disk drops to 
I\ = 0.25. This confirms and quantifies the well-known 
result that low surface brightness disks are sub-maximal 
(de Blok & McGaugh 1997). There is no empirical indi- 
cation that we have reached a lower limit in I\. Disks 
lower in surface brightness than the most extreme consid- 
ered here do exist, though they have yet to be observed 
with sufficient accuracy to be included here. 

The meaning of I\ for high surface brightness galax- 
ies is subject to the caveats discussed in §3.1. In par- 
ticular, the effective definition of maximum disk typi- 
cally accounts for 84% rather than 100% of the veloc- 
ity at the peak of the disk contribution (Sackett 1997). 
If the disks in the highest quartile have 78% of this 
mass, then they contribute 74% of the velocity at 2.2 
scale lengths (M oc V 2 ). This is not very different from 
the sub-maximal contribution of 63% advocated by Bot- 
tema (1993) and is certainly within the galaxy-to-galaxy 
scatter. Kregel, van der Kruit & Freeman (2005) find 
a slightly lower mean velocity contribution, but their 
sample is dominated by intermediate surface brightness 
galaxies, so such a result is consistent with the trend 
apparent in Fig. 7. 

The results in the literature seem broadly consistent, 



bearing in mind that a good deal depends on what is 
really meant by "maximum disk." The most important 
point here is that the degree of maximality of a disk 
depends systematically upon its surface density. There 
is no magic value of that is a fixed fraction for all 
disks (Bottema 1997; de Blok & McGaugh 1996). 

5. CONCLUSIONS 

I have explored the Baryonic Tully-Fisher relation for 
many choices of stellar mass-to-light ratios using a sam- 
ple of high quality data spanning a large dynamic range 
in mass. I provide fits to the BTF for each T*. There 
is a particular choice, based on the minimization of the 
scatter in the local mass-discrepancy — acceleration rela- 
tion (the MDAcc: McGaugh 2004), that also minimizes 
the scatter in the BTF. This optimal BTF is 

M d = 50V/ 

with Vf in km s _1 and the total baryonic mass of a galaxy 
in A4q. This provides a remarkably precise method of 
estimating the mass of rotating galaxies by observation of 
a single global observable, the level at which the rotation 
curve becomes flat. 

The form in which the baryonic mass resides, stars 
or gas, makes no difference to the BTF. Only the sum 
matters. This strongly suggests that the BTF is a funda- 
mental relation between rotation velocity and baryonic 
mass. It further implies that there is no other large reser- 
voir of baryons which matter to the sum: the stars and 
gas observed in spiral galaxies account for essentially all 
of the baryonic mass therein. 

The mass-to-light ratios determined for the optimal 
BTF are in exceptionally good agreement with stellar 
population synthesis models. This consistency, together 
with the agreement between local and global empirical 
relations connecting baryonic mass to the observed dy- 
namics, implies that the baryonic mass is well deter- 
mined. This would appear to solve the long standing 
problem of the uncertainty in the mass of stellar disks. 

Using these robust stellar mass estimates, I have ex- 
amined a variety of evolutionary measures. The stellar 
fraction, mass-to-light ratio, and color all correlate with 
each other as one would expect. Little evolved galaxies 
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with low /* tend to have blue colors and low Y*; more 
evolved galaxies have higher stellar fractions, redder col- 
ors, and higher mass-to-light ratios. These quantities 
are also correlated with disk mass and rotation velocity: 
more massive disks tend to be more evolved. 

The degree to which disks are maximal varies sys- 
tematically with stellar surface density. High surface 
brightness galaxies tend to be more nearly maximal, typ- 
ically with ~ 78% of the maximum disk value at 
[1q = 20.4 mag. arcsec -2 . Low surface brightness galax- 
ies are sub-maximal, with Y* ~ 25% of maximum disk 
at /if = 23.3 mag. arcsec . There is considerable varia- 
tion from galaxy to galaxy. In no case is the stellar mass 



completely negligible at small radii, a fact that is impor- 
tant to mass models and constraints on the inner slope 
of the halo mass distribution (core or cusp). 

The work of SSM is supported in part by NSF grant 
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and Rob Swaters. SSM thanks the Astronomy depart- 
ment of Case Western Reserve University for its hospi- 
tality during a sabbatical visit when much of this work 
was done. 
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